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Subject headings: hydrodynamics — ISM — nuclear reactions, nucleosynthesis, 
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1. INTRODUCTION 

The advent of modern X-ray observatories such as Chandra and XMM-Newton has pro- 
duced a spectacular increase in both the quantity and the quality of the observations of Type 
la supernova remnants (SNRs). Yet, these excellent observations have led only to a modest 
improvement in our knowledge of the physics of Type la supernovae (SNe). Important issues 
such as the nature of the progenitor systems, the last stages of their evolution prior to the 
SN explosion or the physical mechanism behind the explosion itself still remain obscure (see 
Hillebrandt & Niemeyer 2000; Branch & Khokhlov 1995, for reviews). In a recent paper 
(Badenes et al. 2003, henceforth Paper I), we examined the prospects for the identification 
of the explosion mechanism in Type la SNe through the analysis of the X-ray spectra of 
young SNRs. We assembled a grid of Type la SN explosion models, simulated their interac- 
tion with a uniform ambient medium (AM) and calculated the predicted X-ray spectra from 
the ensuing SNRs. The calculated X-ray SNR spectra varied dramatically from model to 
model, demonstrating that it is possible to use young Type la SNRs to probe the details of 
the Type la SN explosion mechanism. 

In this paper, we expand the results that were introduced in Paper I. Our objective 
is to examine the relationship between Type la SN explosions and the X-ray spectra of 
their SNRs within the framework of hydrodynamic, ionization, and spectral simulations. By 
comparing our models with observations, wc aim at improving our understanding of both 
Type la SNe and young, ejecta-dominated SNRs. In § 2, we review the simulation scheme 
used in Paper I, and we discuss the influence of two important parameters which we had not 
hitherto explored: the amount of coUisionless electron heating at the reverse shock and the 
density of the uniform AM that interacts with the ejecta. In § 3 we examine the performance 
of recent 3D Type la SN explosion models in the context of SNRs, and we discuss the ability 
of these 3D models to reproduce the fundamental properties of the X-ray spectra of Type la 
SNRs. Our conclusions are presented in § 4. In order to facilitate the comparison between 
our models and X-ray observations of SNRs, we have generated a library of synthetic spectra. 
This library is presented and discussed in the Appendix. In a forthcoming paper (Badenes et 
al., in preparation), we will make a detailed comparison between our models and the X-ray 
spectrum of the Tycho SNR. 
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2. PARAMETERS AFFECTING THE X-RAY SPECTRUM 

2.1. From SN to SNR: the Simulation Scheme 

Although the X-ray spectra of young Type la SNRs contain much information about 
the structure and composition of the material ejected by the SNe that originated them, this 
information is generally difficult to extract. The ejecta material consists almost entirely of 
heavy elements that are impulsively heated to X-ray emitting temperatures as the reverse 
shock propagates inwards in the reference frame of the expanding ejecta. The propagation 
of the reverse shock is in turn intimately related to the density structure of the ejecta, 
which results in an intricate dynamical behavior of the SNR early in its evolution (see 
Dwarkadas & Chevalier 1998). If there is a significant degree of stratification in the elemental 
composition of the ejecta, different chemical elements are shocked at different evolutionary 
times, after different periods of free expansion, and therefore emit X-rays under different 
physical conditions. This results in a very complex spectrum, which is hard to model and 
interpret. 

The approach taken in Paper I was based on a grid of ID Type la SN explosion models 
that included all the mechanisms currently under debate for the single degenerate Type 
la SN scenario: deflagrations, delayed detonations, pulsating delayed detonations and sub- 
Chandrasekhar explosions. The dynamics of the interaction of each explosion model with 
a uniform ambient medium (AM) of density p^M = 10^^"^ g ■ cm^^ was followed with a ID 
hydrodynamic code. The dynamic evolution of each fluid element in the shocked ejecta 
(i.e., the time evolution of density p and speciflc internal energy per unit mass e), together 
with its chemical composition as determined by the SN explosion model, were used as the 
input to ionization calculations. These calculations included the interactions between ions 
and electrons in the shocked plasma, and they provided time-dependent nonequilibrium 
ionization (NEI) states and electron temperature. Using these ionization states and electron 
temperatures, spatially integrated synthetic X-ray spectra were generated with a spectral 
code by adding the weighted contributions from each fluid element in the shocked ejecta. 
For a more detailed explanation, see Paper I and the references therein. 

Within this simulation scheme, the X-ray spectrum from the shocked ejecta is deter- 
mined by: (1) the density and chemical composition profiles of the SN ejecta from the 
explosion model, (2) the age of the SNR, (3) the amount of collisionless electron heating 
at the reverse shock, and (4) the density of the uniform AM. In Paper I, we analyzed the 
importance of (1) and (2); the impact of (3) and (4) is the focus of this Section. 
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2.2. Collisionless Electron Heating at the Reverse Shock 



The unknown efficiency of collisionless electron heating in SNR shock;s is one of the 
main uncertainties affecting the calculated X-ray spectra of SNRs. Direct application of the 
Ranlcine-Hugoniot relations at the shock front yields 



for each population of particles p, where rup is the particle mass, Vs is the shock velocity and 
k is Boltzmann's constant. Because of the large difference between electron and ion masses, 
the electrons are expected to be much colder than the ions, and the quotient of postshock 
specific internal energies defined as 



is expected to be close to 0, where Zg is the preshock ionization state (i.e., the number of 
free electrons per ion in the unshocked ejecta). However, Cargill & Papadopoulos (1988) 
argued that plasma waves can redistribute energy among cold electrons and hot ions at the 
shock, bringing the value of (5 close to Zg (for a discussion of collisionless electron heating 
see Laming 2000). 

So far, the observational evidence hints at a decreasing level of thermal equilibration with 
increasing shock speeds or Mach numbers in the forward shocks of SNRs (see Rakowski et al. 
2003, and references therein). In the forward shock of Tycho, Ghavamian et al. (2001) found 
a value of Te/Tj < 0.1 by analyzing the optical Balmer emission, while Vink et al. (2003) 
estimated a much lower value at the forward shock of SN 1006 from X-ray observations. The 
only constraint on the value of (3 in the reverse shock of a young Type la SNR comes from 
the absorbed spectrum of the Schweizer and Middleditch star behind SN 1006, where the 
amount of thermal energy deposited in the electrons was found to be negligible (Hamilton 
et al. 1997). The model spectra presented in Paper I were calculated assuming no coUisionless 
electron heating at the reverse shock, effectively setting /3 to the lowest possible value, 
Pmin — Zg ■ me/frii, where is the average ion mass in a fluid element. It is clear from the 
works cited above that, while full thermal equilibration between ions and electrons at the 
shock (i.e., f3 = Zg) is not compatible with the observations, values of (5 larger than (5min 
cannot be excluded. 

The effect of varying amounts of collisionless electron heating at the forward shock on 
the X-ray spectrum emitted by the shocked AM of SNRs in the Sedov stage was discussed in 
Borkowski et al. (2001); here, we shall analyze the impact of a small (but nonzero) amount 
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of coUisionless electron heating at the reverse shock on the properties of the shocked ejecta. 
We will illustrate this effect using a delayed detonation model as an example. Of all the 
classes of one- dimensional Type la explosion models, delayed detonations have been the most 
successful in reproducing the light curves and spectra of Type la SNe (Hofiich & Khokhlov 
1996), and therefore it is of much interest to analyze the details of the X-ray emission that 
these models predict for Type la SNRs. Within the delayed detonations, we chose model 
DDTe because it has the largest amount of intermediate mass elements (Si, S, etc.) in the 
ejecta. This should make it easier to estimate the effects that coUisionless electron heating 
at the reverse shock has on the prominent X-ray lines from these elements. 

Figure 1 shows the shocked ejecta of model DDTe for a uniform AM of density pam = 
10~^^g • cm~^, 430 yr after the explosion (the age of the Tycho SNR). The electron heating 
and plasma ionization processes in the shocked ejecta have been calculated for (5 — Pmin, 
0.01, and 0.1. The unshocked ejecta were assumed to be singly ionized in all cases. As 
discussed in Paper I, the interaction of ejecta with the AM leads to the formation of density 
structures within the shocked ejecta, which strongly affect the distributions of ionization 
states (represented here by the average ion charge Z), electron temperatures Tg, and ioniza- 
tion timescales {t = J riedt). Together with the chemical composition profile of the ejecta, 
these distribution functions determine the spectral properties of each element, and ultimately 
the shape of the emitted X-ray spectrum. An important feature of the shocked ejecta is the 
pronounced density peak towards the contact discontinuity (CD) that appears in all Type 
la SNR models (see Paper I and Dwarkadas & Chevalier 1998). 

In the case with no coUisionless heating (/3 = Pmin), the electron temperature profile rises 
monotonically from the reverse shock to the CD, as internal energy is gradually redistributed 
from the hot ions to the cold electrons through Coulomb collisions. The electron temperature 
profile peaks at the CD, where the fluid elements have been shocked for the longest time and 
have the highest density (the rate at which the ion and electron temperatures equilibrate 
in the shocked ejecta scales with p, see eq. 1 in Paper I). Increasing the value of j3 makes 
the electrons just behind the reverse shock hotter, but the electron temperature drops as 
numerous cold electrons are liberated in the ongoing ionization process and the total internal 
energy in the electrons is redistributed among more particles. For /3 = 0.01, the electron 
temperature proflle eventually relaxes to the proflle obtained without any coUisionless heating 
at the reverse shock, but in the (3 — 0.1 case there is a signiflcant residual temperature excess 
even in the outermost ejecta layers. The average ionization state and ionization timescale 
become severely affected only for /3 > 0.1, when the electrons reach extremely high(~ 10^ 
K) temperatures behind the reverse shock. The ionization process is less efficient at these 
extreme temperatures, leading to lower mean ion charges and ionization timescales in the 
shocked ejecta. 
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The significance of this modification in the electron temperature profile is better un- 
derstood when viewed in the context of the stratification that is inherent to ID Type la 
SN explosion models. In the case of DDTe, as in all delayed detonation models, the inner 
ejecta layers are dominated by Fe and Ni, surrounded by a region rich in intermediate mass 
elements (mostly Si and S, but also Ar, Ca, and others), with O dominating the outermost 
ejecta layers. In Figure 1, this stratified structure has been represented schematically in 
panel b. In this example, the increase in Tg caused by coUisionless electron heating at the 
reverse shock affects primarily the Fe-rich ejecta layers at this age. In terms of the emission 
measure-averaged electron temperature for each element X, (^e)x-i increasing the value of (5 
effectively reverses the approximate ordering in of the ejecta elements that is maintained 
throughout the evolution of the SNR. This is illustrated in Figure 2, which can be contrasted 
with Figure 5 in Paper I. 

In all the calculations presented here and in Paper I, Zg has been set to 1, which is 
generally a good approximation for NEI plasmas in SNRs. Photoionization by UV starlight or 
by X-rays emitted by the shocked material in the SNR could raise Zg, but only by a factor of 
3-4 (see Hamilton & Fesen 1988). After the shock, the heavy element plasma ionizes rapidly, 
and the values of Z{t) in fluid elements that started with different Zg will converge over time. 
Due to this, the postshock X-ray emission is generally insensitive to moderate variations in 
the preshock ionization. In the presence of coUisionless electron heating, increasing Zg for a 
fixed value of (3 raises the number of 'hot' electrons and decreases their temperature in the 
same proportion. X-ray spectra, however, are sensitive mostly to the total internal energy 
transferred to the electrons at the shock and to the final electron temperature, and not so 
much to how the internal energy is distributed among the electrons. For all the values of Zg, 
the flnal electron temperature is similar because the flnal value of Z is very similar. In view 
of this, we do not expect significant deviations from the results presented here for Zg > 1 at 
a fixed value of (5. 



2.3. The Ambient Medium Density 

The density of the AM affects the spectral properties of the elements in the shocked 
ejecta in a dramatic way. This is due to two closely related effects: the acceleration of 
all the collisional plasma processes in denser media on one hand and the scaling of the 
hydrodynamic models with pj^M on the other hand (see section 4.1 in Paper I for a discussion 
of the hydrodynamic scaling). 

These effects are illustrated in Figures 3 and 4, which display the structure of the 
shocked ejecta of model DDTe at an age of 430 yr after the explosion, for a uniform AM 
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of Pam — 5 • 10~^^g • cm~^ and pam — 2 • 10~^^g • cm~^, respectively. These simulations are 
also shown for /? = Pmim 0.01, and 0.1 to facilitate comparison with Figure 1. As a result 
of the scaling laws mentioned in Paper I, a SNR in a denser AM will be in a more evolved 
evolutionary stage at any given time, and vice versa - note how the reverse shock has not 
reached the Fe-dominated region of the ejecta at t = 430 yr for pam = 2 ■ 10"^^g ■ cm~^. The 
mean ionization state in the shocked ejecta, which peaks at Z ~ 10 in the outermost Si- 
dominated layers for pam = 10~^^g • cm~^, rises as high as Z ~ 15 for pam = 5-10~^^g • cm~^ 
in the same region, but only reaches Z ~ 6 for pam — 2 • 10~^^g • cm~^. These differences in 
the mean ionization state correspond to differences in the ionization timescales of roughly 
an order of magnitude throughout the shocked ejecta for each factor 5 increase in pam- The 
electron temperature profiles are affected as well, although only by a factor of ^ 2 — 3. 
The higher densities and faster ionization rates tend to mitigate the effect of coUisionless 
electron heating, favoring the convergence to the canonical /3 = Pmin case (compare panels c 
of Figures 3 and 4) . 

Varying the value of pam has an immediate impact on the emission measures and emis- 
sion measure averaged quantities, mainly through the hydrodynamic scaling laws mentioned 
in Paper I. The approximate scaling of EMx{t) and {T)x{t) is given by eqs. (3)-(5) in Paper 
I: Pa]^^ for the t axis, pam for EMx{t) and p'am fo'^ {'^)x{t)- These approximate scahngs are 
accurate within a factor of 2 for 2 • 10~^^ < pam < 5 • 10~^^g • cm~^, but they might break 
down for values of pam outside this range. The effect of a change of pam on the electron 
temperatures is more complex, and difficult to approximate with sufficient accuracy in view 
of the sensitivity of X-ray spectra to electron temperature. Whereas an accuracy within a 
factor of 2 is reasonable for EMx and tx , which span several orders of magnitude, changes 
by a factor of 2 are too large compared with the more modest (but still up to 2 orders of 
magnitude) range in Tg. 



2.4. Effects on the X-ray Spectrum 

In Figure 5, the temporal evolution of the spectra from the shocked ejecta of model 
DDTe is presented for three values of pam ( 10~^^, 5 • 10~^^, and 2 ■ 10~^^g • cm~^) and two 
values of (3 {I3min and 0.1). A preliminary inspection reveals that variations in pam have 
profound effects on the calculated spectra. As expected from the scaling law for < r >x, the 
plasma ionization state varies greatly, and the presence of different ions results in emission of 
different ionic lines. At pam = 5 • 10~^"^g • cm~^, for instance, the more advanced ionization 
state of Fe leads to a higher flux in the Fe L complex, which blends with O Lya and Mg 
Hea emission at the XMM-Newton CCD spectral resolution. The increase in the Fe Ka line. 
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on the other hand, is due to the higher temperatures in the ejecta. The prominent O Hea 
hne at 0.56 keV, seen at early times for pam = 10~^^g • cm~^, disappears at higher values 
of PaM) because He-like O is ionized more rapidly. The overall higher ionization state of 
the plasma also leads to an increase in the flux of the Lya lines of Si and S, and a shift 
towards higher energies of the Ca Ka line. The shape and flux of the continuum emission 
also change. At lower densities, these effects are reversed. The Fe Ka line and Fe L complex 
virtually disappear, revealing the underlying Ne Hea and Ne Lyct lines at 0.9 and 1.0 keV. 
The O Hea line becomes more important than O Lya, and the Lya and He/3 lines of Si 
and S vanish almost completely, as well as the Ca Ka line. The continuum flattens and the 
emitted flux is generally lower at all energies. 

In contrast with the global effects of variations of Pami changes in the amount of coUi- 
sionless heating at the reverse shock have a different impact on different elements in a model 
with stratifled ejecta like DDTe. For pam — 10~^^g • cm~^, the flux in the Fe Ka line, which 
probes material at higher Tg and lower r than the Fe L complex, is increased by almost two 
orders of magnitude for j3 = 0.1. None of the other elements seems to be affected at this AM 
density, although model DDTe has a significant amount of S, Si and Ca in the inner ejecta. 
This increase in the Fe Ka flux becomes less pronounced with time, and is accompanied by 
a slight change in the shape of the continuum. For pam = 5 ■ 10^^^ g • cm~^, the continuum 
is unaffected and the increase of the Fe Ka line flux is reduced to less than an order of 
magnitude at early times, disappearing completely at late times. At pam — 2 • 10~^^g • cm"^, 
however, the coUisionless electron heating has a more noticeable effect. The shape of the 
spectrum is not changed at low energies, but the flux is somewhat lower at early times for 
f3 = 0.1. At high energies, the level of continuum rises and the flux in the Fe Ka line 
flux greatly increases. The effects of collisional electron heating can be clearly seen even at 
CCD spectral resolution. With a higher spectral resolution such as provided by ASTR0-E2, 
the predicted large temperature increases caused by coUisionless heating (Fig. 2) should be 
detectable through various temperature-sensitive line diagnostics for a number of different 
chemical elements within the shocked ejecta. 

We emphasize that model DDTe is presented here just as an illustrative example. For 
obvious reasons, it is not practical to present the effects of pam and /? on the spectra of all 
the models in our grid. Although the details may vary, the general trends identifled here 
for DDTe can be apphed to most of the other models (for a discussion of other delayed 
detonation models, see Badenes et al. 2005). 

To conclude, we note that coUisionless electron heating at the reverse shock can have 
interesting effects on the spatially resolved X-ray emission. In particular, the enhanced flux 
in the Fe Ka hne discussed above would come mainly from the hotter regions of Fe-rich ejecta 
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close to the reverse shock (see the shape of the electron temperature profile in Figures 1, 3, 
and 4 for values of j3 above 0.01). This scenario is compatible, at least qualitatively, with the 
finding that the Fe Ka emission peaks at a smaller radius than the Fe L and Si Hea emission 
in the X-ray CCD images of both the Tycho (Hwang et al. 1998) and Kepler (Cassam-Chenai 
et al. 2004) SNRs. Collisionless electron heating provides a simpler explanation to the rise 
of the electron temperature profile towards the reverse shock than the relic of an interaction 
with a circumstellar medium invoked by Dwarkadas & Chevaher (1998) for the Tycho SNR. 



3. X-RAY SPECTRUM FROM THREE DIMENSIONAL TYPE la 

EXPLOSION MODELS 

3.1. Type la SN Explosions in 3D: Fiandamental Properties 

In Paper I, we introduced a grid of eight one dimensional Type la SN explosion models 
that included examples of all the paradigms currently under debate: sub-Chandrasekhar 
explosions, deflagrations, delayed detonations, and pulsating delayed detonations. This re- 
duced grid is just a representative sample of a more extensive grid of 19 models, which 
constituted the base for the study of Type la SNRs conducted by Badenes (2004). The re- 
maining 11 grid models are intermediate cases obtained by varying the parameters involved 
in the calculation of each explosion paradigm. We have included these models in the Ap- 
pendix, both for reference in future works and for the convenience of those readers who want 
to use our synthetic SNR spectra for their own research. This grid is one representation of 
our current understanding of one dimensional Type la explosion models, upon which most 
of our knowledge of the physics of Type la SNe is based. 

In view of the recent developments in the field, however, it has become clear that ID 
calculations will soon be superseded by the three dimensional models that have begun to 
appear in the literature (Reinecke et al. 2002; Gamezo et al. 2003; Travagho et al. 2004; 
Garcia-Senz & Bravo 2005). These works have focused on pure deflagrations in 3D, proving 
that they are capable of producing robust explosions, but the ability of these models to 
explain the observations of Type la SNe has not been fully established yet. A common 
feature in all 3D deflagration models, and the most remarkable difference with respect to 
ID models, is the uniform mixing of unburnt C and O material with ^^Ni and the other 
products of nuclear burning throughout the ejecta. This mixing is due to the deformation 
of the flame front caused by Rayleigh-Taylor instabilities, an effect which seems unavoidable 
in 3D deflagrations. There has been some concern that the presence of large amounts of 
C and O in the inner layers of ejecta would lead to a spectral evolution inconsistent with 
optical observations (Gamezo et al. 2003), but complex spectral simulations are required to 
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verify this claim (Baron et al. 2003). Alternatives to the 3D deflagration scenario are being 
explored right now, including delayed detonations in 3D (Garcia-Senz & Bravo 2003; Gamezo 
et al. 2004b, a) and two new explosion paradigms: gravitationally confined detonation (Plewa 
et al. 2004) and pulsating reverse detonation (Bravo et al. 2005), but none of these models 
has been completely understood yet. Here, we study 3D deflagration models in the context 
of the X-ray spectra of young SNRs to provide an independent method of assessing their 
viability for Type la explosions. 

We use a one dimensional average of model B30U, a 3D deflagration from Garcia-Senz 
& Bravo (2005), to illustrate what can be expected from this class of models. The chemical 
composition and density proflle of this model are presented in the Appendix, and they are 
very similar to those of the models obtained by Gamezo et al. (2003) and Travaglio et al. 
(2004), even though the computational techniques and the resolution of the calculations are 
different in each case (see Table 1 in Bravo et al. 2005, for a more detailed comparison of 
these works). This shows that three dimensional deflagrations are relatively well understood, 
and supports our use of model B30U as a representative example of this class. The evolution 
of the emission measures and emission measure-averaged ionization timescales and electron 
temperatures of the principal elements in the ejecta of B30U are shown in Fig. 6 for an inter- 
action with pam = lO^^^g ■ cm^'^ and /3 = /3mm- Interpretation of differences between these 
plots and Figures 4, 5, and 6 of Paper I is not straightforward, because the 3D calculations 
are not fully self-consistent with the ID models of the grid (the effect of the energy deposited 
by the decay of ^^Ni on the density proflle, for instance, has not been taken into account). 
Nevertheless, the main features of the evolution of the shocked ejecta do not depend on such 
details. The most striking property of this model is the similarity in the spectral properties 
((Te) and (r)) of Fe and Si throughout the evolution of the SNR. This is in marked contrast 
to ID models, where the stratiflcation of the ejecta leads to signiflcant differences between 
Fe and Si. The abrupt changes in several plots seen at t ~ 8 ■ 10^° s arc due to the impact 
of the reverse shock on a remnant of unburnt white dwarf material formed in model B30U 
(for details, see Garcia-Senz & Bravo 2005). 

In Figure 7, we show the ejecta spectra of model B30U at the same values of t, p^M 
and (3 as in Figure 5 for model DDTe. The most remarkable properties of these spectra are 
the high Fe L-shell flux and the presence of the prominent Ni Ka line at ~ 7.5 keV (except 
at the lowest AM densities) . This is due to the large amounts of Fe and Ni that are found 
in the outermost layers of B30U, where the density of the shocked ejecta is highest. The 
results are a long Fe ionization timescale, which leads to the enhanced Fe L-shell flux, and 
a high Ni emission measure leading to a strong Ni Ka emission. Such a strong Ni Ka line 
has never been observed in thermal X-ray spectra of SNRs. Another interesting feature is 
the relative weakness of the Si and S Ko; lines. The reason for this is twofold: flrst, 3D 
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deflagrations produce smaller amounts of Si, S and other intermediate mass elements than 
the conventional ID delayed detonations; second, equivalent widths of the Si and S lines are 
smaller because of the strong continuum produced by the large amount of C and O that is 
present throughout the cjecta. These spectral characteristics exhibited by model B30U are 
common to all 3D deflagrations with well mixed ejecta. 

3.2. Compcirison with X-ray observations of SNRs 

We compare the results of our simulations for the 3D deflagration model B30U with 
the basic properties of Type la SNRs. The prediction of similar emission measure- averaged 
electron temperatures and ionization timescales for Si and Fe can be easily tested by exam- 
ining X-ray observations. Wc have searched the literature for young SNRs with published 
good-quality X-ray spectra that have been classified as Type la. Six objects meet these 
requirements: the historical remnants of Tycho, Kepler and SN 1006, and three Large Mag- 
eUanic Cloud SNRs: N103B, 0509-67.5, and DEM L71. We note that the classification of 
Kepler's SNR as Type la is controversial (e.g., Blair 2004). The remnant of SN1006 is not 
suitable for our purposes because it lacks strong Fe emission (Koyama et al. 1995). In the 
case of DEM L71, although the X-ray spectrum of this SNR has been analyzed in some 
detail (see Hughes et al. 2003; van dcr Heyden et al. 2003), wc found no published estimates 
of temperatures and ionization timescales for Fe and Si in the ejecta. The results of the 
analysis of the integrated spectrum for the other four SNRs are summarized in Table 1. 

The spectral properties of all the SNRs considered here show that an important fraction 
of the Fe in the shocked ejecta is emitting under conditions different from those of the Si. 
The authors of the works referenced in Table 1 accounted for this by either adding a spectral 
component made of pure Fc to their fits or by using plane-parallel models that allowed 
to treat Fe and Si separately by assigning different values of Ugt and kT to each element. 
Since the analysis techniques, models, and data quality were different in each case, these 
results can only be compared either with our models or among themselves in a qualitative 
way. Nevertheless, a clear trend can be observed in all four SNRs considered here: the Fe 
component was always hotter than the Si component by at least a factor of 2. The Fe was 
at a lower ionization timescale in three out of four objects: Tycho, Kepler, and N103B. In 
0509-67.5, however, the Fe component has a higher ionization timescale than Si. In this case, 
the statistics of the Fc Ko line were poor, and the ionization timescale of Fc was constrained 
mostly by fitting the Fc L complex. Warren & Hughes (2004) note that their fit to the 
Fe L complex emission was not complete, because a strong line had to be added by hand. 
Improved atomic physics and higher resolution data would be highly desirable to confirm 
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this result for 0509-67.5. 

Since the emission measure averaged ionization timescales and electron temperatures 
of Si and Fe do not differ by more than 30% in model B30U (see Figure 6), we conclude 
that this model is in conflict with the observations hsted on Table 1, at least within the 
limitations of our simulations. This conclusion is extensible to any model in which Fe and Si 
are well mixed throughout the ejecta, and therefore can be applied to all the 3D deflagration 
models for Type la SNe discussed in the previous section. As we have seen in § 2, a plasma 
state with (Te) si < (Te) Fe and (r) si > (t) pe arises naturally in Type la SN models with 
stratified cjccta, such as ID delayed detonations or pulsating delayed detonations, that 
undergo a moderate amount of coUisionless electron heating at the reverse shock. The 
ionization timescales of Fe and Si in 0509-67.5 are clearly incompatible with this scenario, 
but in this case Warren & Hughes (2004) found a very low amount of Fe in the shocked 
ejecta, with Fe to Si abundance ratios below 0.07. While this is very difficult to interpret in 
the context of well-mixed Type la SN ejecta, it could be more easily explained if the reverse 
shock were just entering the Fe-dominated region in stratified ejecta. A detailed comparison 
of our models with this SNR would be required to confirm this hypothesis. 

We emphasize that our simulations based on one dimensional averages are too simple 
to rule out well-mixed 3D Type la SN explosion models. Wc do not account for a number of 
processes that might result in the Fe and Si in the ejecta emitting under different conditions, 
like the Ni bubble effect (Basko 1994; Blondin ct al. 2001) or the formation of clumps in 
the cjccta (Wang & Chevalier 2001). Nevertheless, wc find that the observations of Type 
la SNRs seem easier to explain in the light of Type la SN explosion models with stratified 
ejecta. 

4. DISCUSSION AND CONCLUSIONS 

In this paper, we have examined several important aspects of the X-ray spectral models 
for the ejecta in Type la SNRs that were introduced in Badcncs et al. (2003). We have 
explored the impact of the amount of coUisionless electron heating at the reverse shock, 
(3, and the density of the ambient medium, pam, on the integrated X-ray emission from 
the ejecta in Type la SNR models of different ages. We found that even small amounts of 
coUisionless electron heating can modify the electron temperature profile inside the ejecta 
in a significant way, leading to a region of hot material at low ionization timescales close 
to the reverse shock. In the context of Type la SN explosion models with stratified ejecta, 
this modified temperature profile can affect the emission from the inner layers rich in Fe 
for a broad range of dynamical ages, increasing the fiux in the Fe Ka complex. This could 
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explain why the Fe Ka emission peaks at smaller radii than Fe L in both the Tycho (Hwang 
et al. 1998) and Kepler (Cassam-Chenai et al. 2004) SNRs. The density of the AM also 
has a strong impact on the X-ray emission from the ejecta. For higher values of Pam, the 
SNR is in a more advanced evolutionary stage at a given age, and the ionization timescale of 
the shocked ejecta increases significantly. At lower values of pam, the ionization timescales 
decrease and the evolutionary stage is less advanced. We have provided approximate scaling 
laws to estimate these effects, and discussed their impact on specific emission lines and line 
complexes through an example. 

We have also reviewed the fundamental properties of the recent dcfiagration models for 
Type la SNc calculated in 3D, and their performance in the context of the X-ray spectra of 
SNRs. Using our ID simulation scheme, we have shown that the mixing of fuel and ashes 
throughout the ejecta, which is a common feature of these 3D explosion models, results in all 
the elements in the shocked ejecta of the SNR having very similar spectral characteristics. In 
particular, the emission measure averaged ionization timescales and electron temperatures 
of elements like Fe and Si are always very close to each other. This is in confiict with the 
observations of Type la SNRs in our Galaxy and the Magellanic Clouds, where the Fe and 
Si in the shocked ejecta are found to be emitting under different physical conditions. Within 
the limitations of our ID simulation scheme, these observations are easier to explain with 
Type la explosion models that have stratified ejecta than with models that have well mixed 
ejecta like 3D defiagrations. 

We believe that our models represent a significant improvement over current methods 
of analyzing and interpreting the X-ray emission from the shocked ejecta in SNRs. In order 
to facilitate the comparison between our models and observations, we have built a library of 
synthetic spectra, which is available from the authors upon request. This library is presented 
in the Appendix, where more Type la SN explosion models are introduced, and some aspects 
relevant to the comparison between the synthetic spectra and observations are discussed. A 
detailed example of this kind of comparison in the framework of the ejecta emission from 
the Tycho SNR will be the subject of a forthcoming paper (Badenes et al., in preparation). 

Wc wish to thank Jack Hughes and Jessica Warren for detailed discussions concerning 
0509-67.5. We also acknowledge conversations with Una Hwang and Martin Laming on 
several aspects of the research presented here. Wc arc grateful to the anonymous referee 
for suggestions that helped to improve the quality of this manuscript. This research has 
been partially supported by the CIRIT and MCyT in Spain, through grants AYA2000- 
1785, AYA2001-2360, and AYA2002-04094-C03. CB would like to acknowledge support 
fi-om GENCAT (grant 2000FI 00376) and lEEG in Barcelona, and from grant GO3-4066X 
fi-om SAO at Rutgers. KJB is supported by NASA grant NAG 5-7153. 
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A. A LIBRARY OF SYNTHETIC SPECTRA FOR THE ANALYSIS OF 

EJECTA EMISSION IN TYPE la SNRs 

In this Appendix, we introduce our library of synthetic spectra for the ejecta emission 
in Type la SNRs. The objective of this hbrary is to provide observers with a complete set 
of synthetic spectra for the ejecta emission in SNRs, calculated from an extensive grid of 
Type la SN explosion models, at different values of t, pam, and p. At present, our library 
includes more than 800 synthetic spectra in sequences hke those presented in Figures 5 and 
7 for models DDTe and B30U. For each model, we have generated synthetic spectra for 
several values in the ranges 430 < t < 5000 yr; 2 • 10~^^ < pam < 5 • 10~^^g ■ cm~^; and 
Pmin < /3 < 0.1. In § A.l, we discuss these synthetic spectra in the context of the tools 
that are currently used for the analysis of ejecta emission in SNRs. In§ A. 2 we comment on 
potential applications for our models. In § A. 3, some important caveats and limitations of 
the models are discussed. Finally, in § A. 4, we introduce a number of Type la SN explosion 
models that, together with those presented in Paper I, complete the exploration of the 
parameter space for thermonuclear supernovae. 

A.l. Rationale 

The spectral analysis of the ejecta emission in young SNRs is a complex problem. De- 
spite the spectacular increase in the quality of the observations, it has proved very difficult 
to extract the relevant physical parameters from these observations in a reliable way with the 
available tools. A frequent approach involves the fitting of several more or less sophisticated 
NEI components with varying abundances, electron temperatures, and ionization timescales 
to the observed X-ray spectra (several examples have been cited in § 3.2). The results of 
applying this approach are not unique, and frequently very hard to interpret, because aver- 
age parameters (like Tg or riet) are assigned to a plasma whose physical properties have an 
enormous dynamic range, and where different chemical elements often emit under different 
conditions (see Figures 1, 3, and 4). The determination of elemental abundances in the entire 
volume of shocked plasma, which is crucial for establishing the connection between the SNRs 
and the SN explosions that originated them, is particularly unreliable when it is based on 
this approach. Often, NEI models just provide estimates for the emission measures of the 
elements, Tinder the assTimption of a homogeneous composition, and the difference between 
the ratios of these emission measures and the true abundance ratios in the plasma can be 
several orders of magnitude (see section 4.2 of Paper I). 

The synthetic spectra presented in Paper I and in the present work open new possibilities 
for the interpretation of X-ray observations of Type la SNRs. Without claiming to include 
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all the complex physical processes at play in young SNRs (see § A. 3), these synthetic spectra 
provide a much more accurate representation of the state of the shocked ejecta in young 
Type la SNRs than the simple NEI models currently available within software packages like 
XSPEC. Moreover, since the synthetic spectra are calculated from realistic SN explosion 
models, the connection between the observed spectrum and quantities like the explosion 
energy or the amount of each element present in the ejecta are easy to establish. The 
trade off is that the comparison between our synthetic spectra and X-ray observations is not 
necessarily a straightforward procedure. 



A. 2. Compciring Models and Observations 

Several strategics with varying degrees of sophistication can be followed to compare 
our models to observations. A somewhat crude possibility is to focus on derived quantities 
hke (Te) and (r), as we have done in § 3.2. While this can lead to interesting results, it is 
far better to perform more direct spectral comparisons using the library of synthetic spectra 
that we present here. The most effective way to apply this library will depend on the specific 
observational constraints for the SNR under study. In some cases, like the historical Galactic 
SNRs, the age will be known accurately, but the distance (and hence the total integrated 
X-ray flux and the radius of the forward shock) will be more uncertain. In other cases, 
like the SNRs in the Magellanic Clouds, the distance will be known, but the age will not. 
Reliable independent estimates for pam may or may not be available. In each case, there will 
be more than one way to reduce the dimensionality of the problem. Rather than providing a 
recipe which may not be adequate for a specific case, we make here some suggestions which 
might prove useful in a more general context. 

First, it is important to note that each synthetic spectrum is based on an underlying 
hydrodynamic model, so quantities such as the radius of the forward and the reverse shocks 
and their expansion parameters are available for each spectral model (see Figure 3 and eqns. 
(3), (4), and (5) in Paper I). In principle, it is possible to reverse the problem, find out 
which hydrodynamic models agree better with the observations and thus reduce the number 
of synthetic spectra to consider. In doing so, however, the limitations of ID adiabatic 
hydrodynamics must be considered (see the following section). Second, the selection of 
a particular synthetic spectrum from our library to represent the ejecta component in an 
observed X-ray spectrum may not be trivial. The substantial uncertainties in the atomic 
data and the relative simplicity of the models with respect to real SNRs will probably make 
it impossible to attain a statistically valid fit. Synthetic spectra like ours are more vulnerable 
to these factors, because there is little room for self-adjustment, in contrast to conventional 
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NEI models with variable abundances. If the emissivity of a particular line is underestimated 
in the spectral code used to generate our library, for instance, this cannot be compensated 
by artificially enhancing the abundance of that particular element, as in a conventional NEI 
model. Yet, even if some specific details of the observed spectrum cannot be reproduced, it is 
often possible to find a model whose overall characteristics are in reasonable agreement with 
the observations. Under these circumstances, a procedure needs to be devised in order to 
measure the degree of success of a specific synthetic spectrum. An example shall be provided 
in a forthcoming paper on the Tycho SNR (Badenes et al., in preparation). 

A. 3. Approximations and Caveats 

Our models arc just a simplified representation of the complexity of young SNRs, and 
their limitations have to be considered when making comparison with observations. The 
crucial approximations were reviewed in sections 3.5 and 5 of Paper I, but it is important to 
revisit several issues here. 

The most important simplification is certainly the assumption of spherical symmetry. 
Any description of young SNRs in the framework of 1-D models is necessarily incomplete, 
because it does not include important processes such as ejecta clumping and dynamic in- 
stabihties at the contact discontinuity between shocked ejecta and shocked AM (Chevaher 
et al. 1992; Wang & Chevalier 2001). The degree of ejecta clumping is crucial, and it is 
clear that our ID models (and in particular, the distribution of r and for each element) 
would be invalidated if clumps with a large density contrast like those proposed by Wang 
& Chevalier (2001) were to dominate the emission measure of the shocked ejecta in Type 
la SNRs. In this case, gross inconsistencies are expected to emerge from comparison of 1-D 
models with observations. The degree of ejecta clumping strongly affects the morphology of 
the X-ray emission, and examination of this morphology in Type la SNRs should shed hght 
on this issue. Multi-dimensional hydrodynamical simulations coupled with X-ray emission 
calculations could prove useful for this. 

Another important issue whose impact on the X-ray spectra is hard to estimate is the 
effect of cosmic ray acceleration at the shocks. There is some indication that this process 
might affect the dynamics and X-ray spectra of the shocked AM without significantly mod- 
ifying those of the shocked ejecta (Decourchelle et al. 2000), but more detailed simulations 
are needed to shed hght on this question (see Ellison et al. 2005, for a discussion). 

To conclude this section, we comment on the importance of radiative losses, which have 
received some attention lately in the work of Blinnikov et al. (see Blinnikov et al. 2004, and 
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references therein). In Section 4 of Hamilton & Sarazin (1984), it was shown that radiative 
losses always lead to catastrophic cooling in heavy element plasmas, driving the shocked 
material to infrared- and optically-emitting temperatures. Because no optically emitting 
knots with a composition dominated by heavy elements have been observed in Kepler, Tycho, 
or SN 1006 (see Blair et al. 1991; Smith ct al. 1991), we conclude that radiative losses are not 
dynamically important in young Type la SNRs under usual conditions. Radiative losses are 
not included in our models in a self-consistent way, but we have extended the a posteriori 
monitoring of radiative losses described in section 3.5 of Paper I to the more unfavorable 
case of pam — 5 • 10~^^g • cm~^. Our previous conclusions have been verified: radiative losses 
only affect the outermost layers of the models with the steepest ejecta density profiles. The 
values of trad, as defined in Paper I (the time when the calculated losses exceed 10% of the 
specific internal energy in at least 5% of the ejecta mass) for models with t^ad <5,000 yr are 
provided in Table 2. In these models, our calculations for the properties (density, electron 
temperature, ionization state, and X-ray emission) of the layers that undergo radiative losses 
are not reliable close to or beyond trad, and some amount of infrared or optical emission should 
be expected from this region of the ejecta. The fact that such emission is not observed in 
Kepler, Tycho or SN 1006 suggests that models which predict substantial radiative losses in 
the ejecta are in conflict with observations of these historical SNRs. 



A.4. The Complete Grid of Type la SN Explosion Models 
A.4.i. One Dimensional Models 

Among the eight one dimensional Type la SN explosion models introduced in Paper I, 

one was a sub-Chandrasckhar explosion (SCH), one was a pure detonation (DET), two were 
pure deflagrations (DEFa and DEFf), two were delayed detonations (DDTa and DDTe), and 
two were pulsating delayed detonations (PDDa and PDDe). These explosion paradigms, and 
the details involved in the calculation of the models, are described in Paper 1 (Section 2 and 
Appendix). For the deflagrations, delayed detonations and pulsating delayed detonations, 
the models presented in Paper I were extreme cases, obtained by considering the highest and 
lowest reasonable values of the parameters involved in each calculation. In the case of the 
deflagration models, the relevant parameter is k, which controls the propagation velocity of 
the subsonic flame. For the delayed detonation and pulsating delayed detonation, the pa- 
rameters are ptr, which determines the density at which the transition from deflagration to 
detonation occurs, and l, which determines the flame velocity in the deflagration stage. All 
these parameters are deflned in the Appendix of Paper 1. By varying these parameters, we 
have generated four more deflagrations (DEFb, DEFc, DEFd, and DEFe), four more delayed 
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detonations (DDTb, DDTbb, DDTc, and DDTd), and three more pulsating delayed detona- 
tions (PDDb, PDDc, and PDDd). In Tabic 3 and Figure 8, we present the nucleosynthetic 
output, chemical composition profiles and density profiles of these intermediate models that 
complete the exploration of the parameter space. 



A. 4 -2. Three Dimensional Models 

A brief discussion on the state of the art in 3D calculations of thermonuclear SN explo- 
sions can be found in § 3.1 of this work; for a review see Bravo et al. (2005). Without going 
into the details of how these 3D models are calculated, here we present four one- dimensional 
mappings of 3D models that are representative of the current trends. Their main characteris- 
tics arc given in Figure 8 and Table 4. Model B30U is a 3D deflagration from Garcia-Scnz & 
Bravo (2005), very similar to the models by Gamezo et al. (2003) and Travaglio et al. (2004) 
(see discussion in § 3). Model DDT3DA is a 3D version of the delayed detonation paradigm 
(Garcia-Senz & Bravo 2003). In this model, a detonation was artificially inducted in those 
regions were the fiame resulting from the turbulent defiagration phase was well described 
by a fractal surface of dimension larger than 2.5. We note that this particular model also 
results in very well mixed ejecta, and in fact the properties of its X-ray emission in the SNR 
phase are very similar to those of model B30U. Other delayed detonations in 3D calculated 
with different assumptions for the induction of the detonation result in more stratified ejecta 
(Gamezo et al. 2004b, a). For a comparison between these two kinds of three dimensional 
delayed detonations, see Table 1 and the accompanying text in Bravo et al. (2005). Finally, 
two 3D sub-Chandrasekhar models from Garcia-Senz et al. (1999) have also been included 
in the grid. Model SCH3D0P is a sub-Chandrasekhar explosion calculated in 3D where 
the layer of degenerate He was ignited at one single point, while in SCH3DMP the ignition 
happened at five different points. It is worth noting that none of the 3D models has been 
followed for a sufficient time to account for the effects of the decay of ^^Ni on the density 
profiles. 
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Fig. 1. — Shocked ejecta structure vs. radius for model DDTe interacting with an ambient 
medium of density 10~^^g • cm~^, 430 years after the explosion. The panels show the radial 
distribution of density and specific internal energy (a), mean number of electrons per ion, 
Z, with an indication of the ejecta layers dominated by Fe, Si-S and C-0 (b), electron and 
ion temperatures (c) and ionization timescale (d). The positions of the reverse shock and 
contact discontinuity are outlined by the limits of the temperature plots in panel c. The 
three plots for Z (panel b), (panel c) and r (panel d) represent the values obtained with 
= /?mm (sohd), /? = 0.01 (dashed) and /? = 0.1 (dash-dotted). 
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Fig. 2. — Evolution of (re)x for C, O, Si, S, Ar, Ca, Fe and Ni in the shocked ejecta of 
model DDTe, with p^u — 10"^^ g • cm"^. The top left panel corresponds to (3 — Pmim and 
is the same as panel b of Figure 5 in Paper I, but with a different scale. The top right panel 
corresponds to P — 0.01, and the bottom left panel to P — 0.1. 




Fig. 3. — Shocked ejecta structure vs. radius for model DDTe, with pam = 5-10 ^^g ■ cm 
430 years after the explosion. See Figure 1 for an explanation of the plots and labels. 




Fig. 4. — Shocked ejccta structure vs. radius for model DDTe, with pam = 2-10 g ■ cm 
430 years after the explosion. See Figure 1 for an explanation of the plots and labels. 
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Fig. 5. — Unabsorbed X-ray spectra from the shocked ejecta of the DDTe model, 430, 
1000, 2000 and 5000 yr after the SN explosion, for pam = 10~^'^g ■ cm~^ (top panels), 5 • 
10~^^g ■ cm~^ (middle panels), and 2 ■ 10~^^g ■ cm^^ (bottom panels), convolved with the 
spectral response of the XMM-Newton EPIC MOSl CCD camera. The spectra in the left 
panels correspond to /5 = /3mm, those in the right panels to /3 = 0.1. The Ka lines of Fe, Ca, 
S, and Si, as well as the O Lya line, have been marked for clarity, and fluxes are calculated 
at a fiducial distance of 10 kpc. Note that the spectral code has no atomic data for Ar. See 
the on-line edition for a color version of this Figure. 
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Fig. 6. — EM{t) (top left), < Tg > {t) (top right) and < r > (t) (bottom) for the principal 
elements in the shocked ejecta of model B30U, interacting with a uniform AM of pam = 
10~^'^g ■ cm~^. The crosses in the EM{t) plot represent the total emission measure of the 
shocked ejecta. See the on-line edition for a color version of this Figure. 
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Fig. 7. — Unabsorbed X-ray spectra from the shocked ejecta in the B30U model, 430, 
1000, 2000 and 5000 yr after the SN explosion, for pam = 10~^'*g ■ cm~^ (top panels), 5 • 
lO^^^g ■ cm~^ (middle panels), and 2 • 10~^^g ■ cm~'^ (bottom panels). The spectra in the left 
panels correspond to (3 = Pmin, those in the right panels to /3 = 0.1. See Figure 5 for an 
explanation of the labels and plots. The online edition of the journal contains a color version 
of this Figure. 
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Fig. 8. — Chemical composition and density profiles for the Type la SN explosion models 
that were not presented in Paper I. The abundances represented here are number abundances 
after the decay of all short hfetime isotopes. The density profiles (pnorm) are represented at 
t = 10^ s after the SN explosion, and they have been normalized by pn = 10^^^ g ■ cm~^. See 
the on-line edition for a color version of this Figure. 



Table 1. Spectral properties of Fe and Si in Type la SNRs 



SNR 


Reference 


Age 


Si component ^ 


Fe component 










Spectral model kT log net 


Spectral model 


kT 


log net 






[yr] 


[keV] [cm^s] 




[keV] 


[cm^s] 


Tycho 


Hwang et al. (1998) 


432 


NEI, single Te, 0.86 - 11 

single TLct 


NEI, single Te, 
single net 


> 1.7 


~ 9 


Kcplcr'' 


Kinugasa & Tsuncmi (1999) 


404 


NEI, single Te, 0.77 ±0.06 10.42l[5;'o| 
single net 


NEI, single Te, 
single net 


> 6 


q ^0+0.07 


0509-67.5 


Warren & Hughes (2004) 


< 1, 000 


Plane-parallel 3.13 ± 0.55 9.93 ± 0.02 
NEI shock 


Plane-parallel 
NEI shock d 


10 n+°° 


10.53 ± 0.02 


N103B 


Lewis et al. (2003) 


< 2, 000 


" Plane-parallel ^1 > 12 
NEI shock 


Plane-parallel 
NEI shock 


> 2 


~ 10.8 



'^In all the cases listed here, the 'Si component' also included some Fe, as well as other elements. 
''The identification of this SNR as Type la is controversial, see Blair (2004). 

'^Kinugasa & Tsunemi (1999) also fitted the spectrum with an NEI component for the Fe Ka line plus a more sophisticated 
model for the rest of the shocked ejecta and AM, based on the self-similar solutions of Chevalier (1982), coupled to NEI 
calculations in a plasma with homogeneous abundances. However, they do not give explicit values for kT and net of the reverse 
shock in this model. For simplicity, we use the results of their two-component NEI model, which gives a fit of similar quality. 

'^We list the results for the best fit model in Warren & Hughes (2004), which includes a nonthermal continuum (their model S). 
Assuming a thermal continuum (their model H), these authors obtain similar results; kTsi = 2.23±0.29, lognetsi = 9.94±0.02, 
kTpe = 10.00t^i4, lognetpe = 10.53 ± 0.02. 

"Hughes et al. (1995). 
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Table 2. Values of trad 



Model pam 


= 10 ■ cm 


Pam = 5-10 ^^g ■ cm ^ 


DEFa 


3.0 • 10^° s 


2 • 10^° s 


DEFc 


2.7-10^° s 


1.6 ■ 10^° s 


DEFf 


2.4 • 10^° s 


1.2 ■ 10^° s 


PDDe 




1.5 • 10^^ s 



'^Only models with trad <5,000 yr (1.58 ■ 10^^ s) are listed. The values of trad for other DEF 
models (DEFb, DEFd and DEFe) are comparable. All calculations were done with /3 = Prnin 
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Table 3. Properties of the additional Type la ID explosion models 



Model 


Pcira- 


ptr 


Mqj ecta 




■f^raax ^ 


A»15 b 




Mc+o 


Msi 


Ms 


MAr 






meter ^ 


[g • cm~^] 


[M©] 


[10^1 erg] 


[mag] 


[mag] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


DEFb 


0.08 




1.37 


0.64 


-19.14 


0.94 


0.61 


0.61 


0.025 


0.017 


0.0040 


0.0043 


DEFc 


0.10 




1.37 


0.74 


-19.29 


0.99 


0.68 


0.55 


0.021 


0.014 


0.0032 


0.0032 


DEFd 


0.12 




1.37 


0.80 


-19.34 


1.02 


0.71 


0.52 


0.021 


0.014 


0.0032 


0.0034 


DEFe 


0.14 




1.37 


0.81 


-19.29 


0.98 


0.73 


0.49 


0.021 


0.013 


0.0029 


0.0028 


DDTb 


0.03 


2.6 ■ 10^ 


1.37 


1.36 


-19.67 


1.11 


0.98 


0.05 


0.10 


0.084 


0.022 


0.027 


DDTbb 


0.01 


2.5 ■ lO'^ 


1.37 


1.31 


-19.66 


1.12 


0.99 


0.05 


0.10 


0.084 


0.022 


0.027 


DDTc 


0.03 


2.2 ■ lO'^ 


1.37 


1.16 


-19.51 


1.11 


0.80 


0.12 


0.17 


0.13 


0.033 


0.038 


DDTd 


0.03 


1.5 • 10^ 


1.37 


1.08 


-19.30 


0.94 


0.72 


0.14 


0.20 


0.15 


0.037 


0.043 


PDDb 


0.03 


2.2 • 10'^ 


1.37 


1.36 


-19.72 


1.14 


1.04 


0.03 


0.085 


0.070 


0.018 


0.022 


PDDc 


0.03 


1.5 ■ 10^ 


1.37 


1.25 


-19.64 


1.11 


0.98 


0.04 


0.11 


0.093 


0.024 


0.029 


PDDd 


0.03 


1.2 • lO'' 


1.37 


1.24 


-19.53 


1.04 


0.89 


0.05 


0.15 


0.13 


0.034 


0.041 



^The parameter given is k for the DEF models and t for the DDT and PDD models (sec the Appendix of Paper I for details). 
''The values of Mmax and for the light curves were calculated by I. Dominguez (private communication, 2003). 
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Table 4. Properties of the 3D Type la explosion models 



Model 






Mfo 


Mc+o 


Msi 


Ms 


MAr 


Mca 




[Mo] 


[10^1 erg] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[M©] 


[Mq] 


B30U 


1.37 


0.42 


0.53 


0.66 


0.045 


0.011 


0.0019 


0.0017 


DDT3Da 


1.37 


0.78 


0.76 


0.38 


0.063 


0.027 


0.0066 


0.0072 


SCH3DOP 


1.02 


1.14 


0.58 


0.23 


0.064 


0.035 


0.0093 


0.0077 


SCH3DMP 


1.02 


1.19 


0.67 


0.07 


0.081 


0.054 


0.019 


0.017 



